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erator for the benchmark process e^e~ — > hadrons. Results are presented at parton and 
hadron levels, demonstrating a smooth merging between the matrix elements and parton 
shower and an improved description of LEP data. 
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1. Introduction 

Monte Carlo event generators have been successful in providing a full simulation of the 
physics at collider experiments and have proven to be invaluable tools in both planning 
future experiments and analysing data from current experiments. They use a flexible 
and convenient event-by-event description which allows a range of physics models to be 
implemented and as such they describe a wide variety of phenomena. In particular, they 
provide a means, via the parton shower, of evolving from hard scales, where partons are 
produced in fixed-order perturbation theory, to soft scales where non-perturbative models 
must be applied. 

Besides the simple fixed-order matrix elements which describe the hard processes, 
the key perturbative element of Monte Carlo event generators is the parton shower. In 
this phase of the simulation, partons undergo DGLAP evolution from some initial scale. 
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characteristic of the underlying hard scattering, down to the hadronization scale. In Her- 
wig++ [1,2] the evolution variable of the parton shower is chosen such that branchings are 
ordered in their opening angles, accounting for the effects of QCD coherence [3]. This cor- 
responds to resumming next-to-leading logarithm (NLL) terms to all orders in perturbation 
theory. 

While the parton shower accurately simulates soft and collinear radiation, it does not 
provide a reliable description of hard (high transverse momentum) emissions. In fact, 
requiring that the regions of phase space into which each shower progenitor (the partons 
which initiate the shower) can emit do not overlap, to avoid double counting, gives rise to 
a dead zone: a region of phase space corresponding to high transverse momentum, wide 
angle, emissions which none of the showering particles can emit into. Even within the 
accessible shower regions of phase space, the distribution of radiation involves some degree 
of approximation, since at any given fixed order in perturbation theory, the parton shower 
effectively approximates the real emission corrections to the hard scattering process by a 
product of splitting functions and Sudakov form factors, summed over all combinations of 
branchings which give rise to the same final state. These approximations account for the 
NLL corrections associated with soft and collinear radiation in the perturbative scries. 

In order to improve the description in the large transverse momentum region, the 
parton shower can be combined with exact fixed-order matrix elements. The earliest and 
simplest means of forming this combination is known as the matrix element correction 
method [4]. This corrects the hardest emission generated by the parton shower such that 
it is distributed according to the real single emission matrix element squared. This tech- 
nique has been successfully applied to important processes in a number of generators [5,6], 
including Herwig+-|- [2]. 

In recent years more general matrix element-parton shower merging algorithms have 
been introduced. These combine tree-level matrix elements with parton showers, for a 
given process, for all parton multiplicities below some maximum N. Hence these algo- 
rithms correct all distributions involving up to N external partons, instead of just that of 
the hardest emission. Several schemes of this type have been developed and successfully 
implemented in event generators. The most well known of these are the CKKW [7-9], 
CKKW-L [10], MLM [11] and pseudo-shower [12] methods. All these methods have the 
same general approach [13, 14] whereby the phase space for parton emissions is divided into 
two regions by a merging scale yj^g, defined in some jet measure. Above the merging scale, 
emissions are described by exact matrix elements while below it emissions are produced by 
the parton shower. 

In this article we present a matrix element merging scheme based on the CKKW 
algorithm. A fundamental ingredient in the CKKW method is the association of a pseudo- 
shower history to the configurations generated according to the fixed-order matrix elements. 
Each shower history is constructed by clustering the two most closely separated partons, 
according to the transverse momentum measure defining the merging scale, until a leading- 
order parton configuration is obtained. The resulting branchings in the shower history are 
therefore ordered according to the jet measure, which may not equate to the ordering 
variable of the parton shower, as is the case for the angular-ordered parton shower of 
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Herwig++. This discrepancy is understood to give rise to serious problems, in particular 
it spoils the colour coherence properties of these shower algorithms. Although this was 
already noted, and an attempt made to address it, in the original CKKW paper, realisations 
of the method highlight the fact that the colour structure in the events is nevertheless in 
conflict with that expected on the grounds of colour coherence, moreover, they show that 
this is not simply an esoteric consideration but a cause of significant practical problems, 
including a dependence on the unphysical merging scale [12, 13]. 

We shall present and validate a modified version of the CKKW method, intended to 
optimise the implementation of these colour coherence effects, by a fully consistent merging 
of an angular-ordered parton shower with fixed-order matrix elements. The idea behind 
this method was originally proposed by Nason in Ref. [15]. The central result of that 
theoretical work is the observation that the parton shower may be formally decomposed in 
terms of truncated showers, hard emissions, and vetoed showers. Reference [15] advocates 
that the CKKW algorithm may then best model the coherent emission of radiation by 
including these truncated showers, consisting of only soft emissions, prior to, and between, 
the hard emissions in the shower history, thereby rendering it angular-ordered. In the 
following we will develop the full details necessary for our practical implementation of this 
idea for the process e~^e~ — > hadrons and compare the results of it to LEP data. 

This work builds on the infra-structure introduced in Refs. [18, 19] for the Herwig-|— |- 
angular-ordered parton shower, where the truncated shower was initially developed for 
use in the POWHEG [15-17] next-to-leading order matching scheme applied to Drell-Yan 
vector boson production. This was the first time a full implementation of the truncated 
shower had been developed^. 

While this work was in production, Ref. [20] appeared on the arXiv, employing sim- 
ilar techniques to those described in the present work for matrix element merging with a 
transverse-momentum-ordered dipole shower. 

The paper is set out as follows. In Sect. Q we review the original CKKW merging 
prescription. In Sect. ^ we go on to describe the way in which the angular-ordered parton 
shower may be decomposed into hard emissions, truncated showers, and vetoed showers. 
Having introduced the relevant conceptual ingredients we then give a more detailed tech- 
nical description of our modified CKKW algorithm in Sect. ^. In Sect. ^ we present a 
validation of our algorithm by comparing to LEP data for e~^e~ — > hadrons, before giving 
our conclusions in Sect. |6[ 

2. CKKW merging 

In this section we present an overview and discussion of the original CKKW algorithm 
for the process e~^e~ hadrons. We first describe the algorithm for the case where the 
parton shower evolution variable is identical to the merging variable before describing the 
adaptations which must be made for the Herwig-|-+ angular-ordered parton shower. 

^An approximate, single emission truncated shower is also described in Refs. [21,22]. 
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2.1 Transverse- momentum-ordered CKKW merging 

The algorithm is simplest if the merging variable is the same as the ordering variable of 
the parton shower. We therefore consider first the case where we have a single transverse 
momentum variable q as the parton shower evolution and merging variables. 

The basic principle underlying the CKKW approach, is that the distribution of ra- 
diation in the region of phase space where all partons are separated by an amount q 
greater than the merging scale qMS^ should be given by tree-level matrix elements, while 
for q < qMS it should be given by the parton shower. The algorithm then requires, as 
input, samples of events of the process with up to N partons in the final state. These 
input samples are easily obtained using fully automated tree-level event generators such as 
Madgraph/MadEvent [25,26]^. As well as producing the events themselves, for each sample 
with n partons the generator will provide a finite, tree- level, jet cross section a^^^ (Qms)- 

Naively, with the input events in hand, one might then consider filling the remaining 
phase space by selecting events from each sample with n partons, with a probability pro- 
portional to a^^\ and simply invoking the parton shower on each of the external legs, 
starting from the scale qj^.jg- However, the merging scale g^^^ is not a physical parameter 
and so all distributions of partons should be insensitive to its value. This would certainly 
not be the case for such a naive procedure, since the distribution of radiation from the 
parton shower and the fixed-order matrix elements are known to differ, especially in the 
regions corresponding to high and low q emissions. The great success of the CKKW al- 
gorithm is in its ability to correct for the mismatch at the phase space partition g^^^ by 
providing a smooth, physical, interpolation between the matrix element distribution at 
high q values and that of the parton shower in the low q region. 

To illustrate how this works, consider the simplified case of merging only samples 
of 2- and 3-parton events, with q > q,^jg, for e~^e~ hadrons, with a q'-ordered parton 
shower. In general, the parton shower cross section analogous to o"!*^^^ (Qms)^ with n 
partons resolved at the merging scale, may be written as the product of the leading-order 
cross section together with a set of Sudakov form factors and splitting functions. The 
product of these splitting functions and the leading-order cross section approximate the 
exact tree-level n-jet cross section. For the case of three partons this cross section is 

(PS) 2 f^' 

(.(1ms) = (^2x2 [Ag / dqas (g) T^-^gg (q) Ag(g,g^,s), (2.1) 

-^Ims 

where qj is the scale at which the parton shower is initiated and OisiQ)^fj^ijiQ) is the 
probability for a parent parton ij to branch into two daughter partons i and j, in the interval 
[q,q -(- dq] ^. The overall normalisation factor £72 is simply the leading-order cross section. 
Finally, in Eq. (^), Ar^.(g, 5^,^^) is the Sudakov form factor, which can be interpreted as 
the probability for the parent parton ij to evolve from a scale q down to the scale g^^^ 

^Currently, computational efficiency limits the total number of final-state particles to around six. 
^The dependence on auxiliary splitting variables has been suppressed. 
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without undergoing a resolvable branching, 



^riiQ,QMs) =exp 



^ r d,'a.(,')r~.^,,(<z') 



(2.2) 



In parton shower (NLL) expansions of the jet cross sections, such as that in Eq. (p.lj), 
the exact tree-level matrix elements are approximated by the product of the leading-order 
cross section and splitting functions. In order to improve the parton shower with exact 
tree-level matrix elements, this product should be replaced by the corresponding exact, 
tree-level jet cross section. 

The CKKW merging should not affect the NLL expansion of the jet cross section 
therefore the NLL expansion of the matrix element contribution should give the result in 



Eq. (2.1). Since a NLL expansion of the tree-level matrix elements yields a corresponding 
product of parton shower splitting functions, it is clear that in order to retain the NLL form 
of Eq. (^]^), the matrix element contribution above Qj^jg should be given by configurations 
generated according to the tree- level jet cross sections reweighted by appropriate Sudakov 
and running as factors. 

In order to determine appropriate reweighting factors for events from the tree-level 
generator, a pseudo-shower history must be assigned to each event. This shower history 
interprets the set of external parton momenta as a set of branchings originating from a 
leading-order configuration. This procedure gives rise to a set of nodal values, qi, for the 
scales at which each pseudo-branching occurred. These scales provide the arguments for 
the Sudakov form factors and as factors with which the configuration should be reweighted. 
In the original CKKW publication, this pseudo-shower history is assigned by repeatedly 
clustering the pair of partons^ with the smallest separation according to the jet resolution 
variable, until only the particles of the leading-order process remain. 

In the case being considered, where the evolution variable has been taken to match the 
merging scale variable, combining the matrix elements with the parton shower is straight- 
forward. The parton shower evolution can be split into two parts: first an evolution from 
the initial scale down to the merging scale q]^Jg', then an evolution from the merging scale 
down to the hadronization scale qq. This results in a simple procedure for attaching the 
parton shower to the reweighted matrix elements, where each external parton produces a 
shower line evolving from the merging scale. 

The full CKKW algorithm then proceeds as follows: 

1. a jet multiplicity n is generated with probability 

(ME) , X 

p _ \1ms) /n q^ 

where all cross sections are evaluated at a fixed strong coupling asj^j^] 

2. a configuration of n parton momenta is generated according to da^^^\qj^^g); 



*Only pairs of partons whose flavours correspond to allowed branchings are considered. 
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3. external partons are clustered, defining a pseudo-shower history with a set of nodal 
scales qi\ 



4. the configuration is reweighted by Sudakov and as factors: each internal line between 
two nodes at qi and gj+i contributes a factor of Aj(gj, g^^c,)/Aj(gj+i, 5^,^^), each 
external line emanating from a node with scale qi contributes A j (gj , qj^^^ ) , while each 
node itself contributes "^^'^'^ ; 

5. the parton shower is invoked on each external parton from a starting scale of qj^jg- 

This scheme is independent of the merging scale to NLL order [7] . We have reweighted 
configurations such that the NLL three-jet cross section resolved at the merging scale is 
given by Eq. ( p.ll ). This NLL cancellation of merging scale dependence can be seen by 
considering the cross section for three jets resolved at the hadronization scale. This cross 
section is given by the sum of the probability of generating a single emission in the matrix 
element region and none in the parton shower, together with the probability of generating 
no emissions in the matrix element region and a single emission in the parton shower region. 
The cross section is 

^+^^^) (go) = 4*'''^ ilMs) [A, (g.«, go)]' A, (g,,„ go) (2.4) 

+ a2 X 2[Ag{qi,qo)] / dqas{q)Tg^gg{q) Ag{q,qo), 
J go 

where a^^^^ (Ims) reweighted matrix element contribution for three jets resolved 



at the merging scale. The first term in Eq. (^^) corresponds to a single emission above 
the merging scale followed by parton shower evolution from the merging scale down to 
the hadronization scale with no resolvable emissions. The second term corresponds to no 
emissions above the merging scale followed by a single parton shower emission below the 
merging scale. In the NLL expansion of Eq. ( |2.4| ), we replace a'^^^ (Qms) NLL 



parton shower approximation in Eq. (pTf). This results in a simplification of Eq. (2.4) 



2 

0-3 (go) = 0-2 X 2 [A, (g/, go)] / dqas{q)Tq^ggiq) Ag{q,qQ), (2.5) 

■J go 

yielding the expected NLL parton shower cross section for a single resolved emission which 
is independent of the merging scale. 

2.2 Angular-ordered CKKW merging 

The merging variable used to define the jet cross sections must regulate both soft and 
collinear singularities, so it must be a transverse momentum measure. The merging variable 
in the original CKKW publication is defined in terms of the Durham jet measure [23] for 
two partons i and j, 

2min (Ef,E^) 

= ^ (1 - cos 9ij) , (2.6) 
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where Eij are the energies of the two partons, Otj is the angle between the two partons 
and s is the centre-of-mass-energy squared. The merging transverse momentum variable is 
defined by 

k± = ^/yijS. (2.7) 

The parton shower with which we wish to merge the matrix elements may not be 
ordered in transverse momentum, in which case the merging variable cannot be chosen to 



be the same as the evolution variable, as was assumed in Sect. 2.1 



In the Herwig++ parton shower, splittings are described by the variables {q,z,(j)), 
where q is an angular-ordered evolution variable, z is the momentum fraction of the emitted 
parton and (j) is the azimuthal angle of the branching. The evolution variable is defined by 



z^(l - z ff =p]_ + {l- z)m^ + zrrij - z{l 



z)m~., 



(2.1 



where p_L is the relative transverse momentum of the branching and rar-. and mjj- are the 
masses of the parent and child partons respectively. 

The probability for a branching ij — > ij in the phase space measure q — > g + dq" is given 

by 

i2 



dP~. ..{q,z) 



27r 



— TT^zP- ■■(z,q), 
q2 



where Pn 



) is the corresponding quasi-collinear splitting function [27]. The Sudakov 
form factor, giving the probability of a parton shower line of flavour ij evolving from a 
scale qi down to q2 undergoing no resolvable emissions, is given by 



exp 



92 



dV^- ■■(q,z) 



(2.10) 



In order to accommodate the fact that the evolution and merging variables are not 
identical, the CKKW algorithm must include some additional features to that outlined 



in Sect. 2A. Changes must be made to the Sudakov form factors with which the matrix 
elements are reweighted and the initial conditions with which the parton shower is invoked, 
furthermore, when the shower is invoked, a veto must be applied to prevent it generating 
emissions with k_i_{q,z) > k±j^j^. 

The Sudakov form factor used for the matrix element reweighting, defined in Eq. (^^) , 
corresponds to the probability of evolving from a scale q down to the hadronization scale 



with no emissions resolvable at the merging scale. In the case of Sect. 2.1, this was achieved 



by setting the lower limit on the integral to q,^,ig, however, now this cut, defining what is 
meant by a resolvable emission, must be implemented as a ^-function in the Sudakov form 
factors used in step ^. The Sudakov form factors for the reweighting are then given by 



^p;kA 



exp 



~. . . 92 



(2.11) 



The prescription for constructing the Sudakov weights is then identical to that in Sect. 2A 
except for factors of z in the scale from which each child evolves, which are required for the 
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angular-ordered evolution. Each intermediate line, connecting branchings at {qi,zi) and 
{Q2,Z2) in the pseudo-shower history, contributes a factor 



Each external line, from a branching at (q, z) in the pseudo-shower history, contributes a 
factor 



2.3 Highest multiplicity treatment 

The original CKKW publication did not treat the highest multiplicity matrix element 
contribution any differently to the other multiplicities. In Refs. [9, 12] it was noted that 
a different treatment of highest multiplicities must be employed in order to fill the phase 
space in the matrix element region to all orders in as- Since computational limits mean 
that only matrix elements with up to a maximum of N partons can be calculated, the 
standard approach leads to a maximum of N partons being generated above the merging 
scale. The parton shower generates to all orders in as and therefore we should also let the 
matrix element region generate to all orders. This can be achieved by allowing the highest 
multiplicity channel parton shower to generate emissions in the region with k± less than 
that of the lowest transverse momentum of the matrix element emissions, k±j^. This is 
achieved by changing the scale of the parton shower vetoes and Sudakov form factor cuts 
from k±j^,jg to k±j^ . 

2.4 Problems with the algorithm 

The above procedure is heavily reliant on having an exact mapping between the shower 
variables and the merging measure k±j^jg so that the parton shower vetoes and Sudakov 
cuts can be correctly applied. A mapping from the momentum clustered in step ^ to the 
corresponding shower variables is also required, so that the correct scales for the Sudakov 
reweighting and initial shower conditions are obtained. In practice obtaining such mappings 
may be difficult due to the complexity of the shower kinematics. 

The initial scale at which the parton shower is invoked is vital to the algorithm. Ini- 
tiating the parton shower directly from the merging scale would result in a radiation gap, 
where emissions with transverse momentum less than the merging scale but evolution scale 
greater than the merging scale are missed. In the angular-ordered shower, this radiation 
corresponds to soft, wide-angle emissions. The original CKKW publication attempts to re- 
solve this by invoking the parton shower from each external parton at a scale corresponding 
to the node at which it was 'created' in the pseudo-shower history. Although adopting this 
maximal initial scale helps fill the radiation gap, the extra soft, wide-angle radiation that 
results, is emitted from the external parton in the pseudo-shower history, rather than the 
intermediates, as implied by colour coherence [15]. The original CKKW publication argued 
that this should be a sub-leading effect, however, it will certainly change the colour struc- 
ture of the configuration, which may cause problems when non-perturbative hadronization 
models are applied. 






(2.13) 
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The original CKKW algorithm also assumes that the clustering of momentum in step 
m and subsequent mapping to parton shower variables results in a set of emission scales 
that respect the ordering of the parton shower, i.e. 

qi> qi> ... > qn> %■ (2.14) 

The clustering scheme of the original algorithm does not guarantee this. 

These problems were studied in Ref. [13] and found to result in discontinuities around 
the merging scale and a false dependence on the merging scale when the parton shower was 
not ordered in transverse momentum. In Ref. [12], a study of the algorithm with angular 
and virtuality-ordered parton showers was presented. In that work, a number of ad hoc 
adaptations were applied and tuned in order to achieve a reasonably smooth merging at 
the parton level, nevertheless, some problems remained at the hadron level. In this article 
we aim to overcome these problems with a set of well motivated modifications based on 
the POWHEG shower reorganisation. 

3. Shower reorganisation 

In the POWHEG next-to-leading-order matching scheme it was shown that a general par- 
ton shower may be rearranged such that the hardest emission^ can be generated first. The 
hardest emission can then be generated according to the exact NLO matrix elements such 
that inclusive observables are distributed according to the NLO cross section while the 
soft /collinear resummation of the shower is undisturbed. In the POWHEG shower reor- 
ganisation, the hardest emission is then dressed by inserting a truncated shower of wide 
angle, soft emissions prior to it, and a vetoed shower consisting of smaller p±, smaller angle 
emissions, after it. 

The CKKW algorithm generates a set of n emissions above the merging scale y^^jg 
according to exact tree- level matrix elements up to 0{ag). This defines a set of n hard 
emissions. In order to reproduce the full shower around this set of hard emissions we employ 
a generalisation of the POWHEG shower reorganisation. In the following a review of the 
POWHEG reorganisation is presented followed by its extension to the CKKW case. The 
notation used relates specifically to that of the Herwig-|— |- shower, however the treatment 
is independent of the details of the parton shower. 

3.1 POWHEG reorganisation 

A general Herwig-|— |- parton shower, evolving from a parton of flavour ij and initial scale 
qi, can be represented by a generating functional S^-[qi). The evolution of the parton 
shower may be expressed by the recursive equation, 

S.^{qi) = A.^iqj, qo)S..iqo) + f ' A.^{qi, q)dV.^^^^{q, z)S, (zq) S, ((1 - z)q) . (3.1) 

The first term in Eq. ( p.lj ) is a no emission term corresponding to the probability of 
evolving from the initial scale to hadronization scale with no resolvable emissions generated. 

^The hardest emission here refers to the emission with highest transverse momentum. 
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The second term represents the probabihty of producing at least one emission, with the 
first generated at a scale q and further showers evolving down to the hadronization scale. 

It is possible to expand Eq. (^) to explicitly show the hardest emission of the shower. 
The hardest emission is described by the shower variables {(jh, ^h, <A/i) and has an associated 
transverse momentum p_L,^ . The shower may produce any number of other emissions before 
the hardest emission and any number of emissions after it but all of these must have 
P± < P±h ■ '^^^ shower can therefore be written as 



11 



10 



S-j {qi,qh;p±J<iVT,_ 



{qh,Zh)Sl izhqh;P±^) SJ ((1 - Zh)qh;p±^ 

(3.2) 



where S'^ refers to a truncated shower and refers to a vetoed shower. The truncated 
shower in Eq. (^^) is responsible for evolving from the initial scale down to the scale of 
the hardest emission producing any number of emissions with transverse momentum less 
than p±^. The emissions within the truncated shower correspond to soft, wide angle gluon 
emissions, and so do not change the flavour of the shower line. The vetoed shower evolves 
from the scale of the hardest emission down to the hadronization scale also generating only 
emissions with transverse momentum less than its evolution is defined by 



Sl {qh;p±J = ^~j{qh, qo)S~jiqo)+ / \-,{qh, g)dP^._, .(g, z)0{p 

Jqo 



±h-P±iq,z)) (3.3) 



xSY {zq]p±^) Sj ((1 - z)q;p±J . 



The recursive equation describing the evolution of the truncated shower is given by 



S.-. {qi,qh;p±J = ^ij{qi,qh)S-;j{qh)+ 



Ql 



^dqi, q)dVr.^,^iq, z)e{p^^ - p^iq, z)) (3.4) 



xS~. {zq,qh;p±i,) {{1- z)q;p_i_J . 



The Sudakov form factors and splitting functions appearing in the Eqs. (3.4, 3.3) 
are identical to those in the standard shower equation of Eq. ( |3.1D with the exception 
that the splitting functions in both new showers have an additional ^-function. This 9- 
function guarantees that no emissions with transverse momentum greater than that of 
the hardest emission are generated. Standard Monte Carlo techniques require that the 
splitting functions of a parton shower match those appearing in the Sudakov form factors. 
The introduction of the ^-functions mean that this is not the case for the vetoed and 



truncated showers in Eq. (3.2), we highlight this in our notation with a bar. 

In order to make the truncated and vetoed showers suitable for a Monte Carlo treat- 
ment, the original POWHEG publication [15] splits the Sudakov form factor appearing in 
Eqs. (p.4|, |3.3|) into two parts according to 



A 



Zit 



,qi+ij 



Aj {ziqi , qi+i ; p^JAf{ziqi, qi+i ; p±^ ] 



(3.5) 
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Here, Aj refers to a vetoed Sudakov in which the exponent contains a 0-function, which 
matches that in the sphtting function of Eqs.(^!J, |3.3| ) and is given by 



'^fAziqi,qi+i;p±J =exp 



E / . {q, z)e {p^, - ps_ {q, z)) 



(3.6) 



The other factor, Aj , contains the opposite 0-function and is referred to as a remnant 
Sudakov given by 



^fj{ziqi,qi+i\P±^) =exp 



91+1 ~. 



(3.7) 



The combination of the sphtting functions in Eqs. (3^, 3^) and the vetoed Sudakov form 
factors result in a parton shower that may be generated with standard vetoes allowing 
only emissions with px_ < p±f^, however, the presence of the remnant Sudakov form factors 
appears to spoil this picture. On the contrary, it turns out that the seemingly awkward 
remnant factors have a key role to play in formalising how to generate the hardest emission 
first. In Ref. [15] it was proven that the hardest emission in the shower is generated 
along the hardest line, the line for which all Zi > and, moreover, all non-soft emissions 
preceding it give rise to subleading collinear contributions. This means that the truncated 
shower may be considered as comprising solely of soft gluon emissions and that Zi can be 
effectively replaced by one in all of the associated remnant factors. It is also shown that 
the 0-function in the remnant Sudakov form factor exponent is zero for scales less than q^ 
and so the replacement Zi ^ 1 also holds for the vetoed showers. The net result of these 
replacements is that the product of all remnant Sudakov form factors combine as a single 
remnant Sudakov factor: 



A|(Q'/,go;p±, 



exp 



91 



90 



{q, z)9 {p±{q,z) -p±,J 



(3.8) 



The POWHEG treatment which has been outlined, results in a reorganisation of the 
shower such that the hardest emission may be generated first. The Monte Carlo interpre- 
tation of this reorganisation is: 



1. the hardest emission {qh,Zh,(j)h) is generated according to the appropriate splitting 
function reweighted with the remnant Sudakov factor of Eq. (|3.8| ) ; 



2. a truncated shower, allowing only non-flavour-changing emissions with p± < p±^ is 
initiated, evolving the shower from qj down to ijh', 

3. the hardest emission is forced with shower variables {qh, Zh, (ph)', 



4. showers with a veto, allowing only emissions with p± < p±f^, evolve all external lines 
down to the hadronization scale. 
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A shower generated in this way should differ from the standard shower by only sub-leading 
terms. 

The vetoed shower is defined by 

rqh 

S}Aqh;V±n) = ^}A^hAo)S}m+ / Ag(g,,g)dP~.^.^.(g,z)e(p^,-p^(g,z)) (3.9) 

xSY izq;p±J 5]^ ((1 - z)q;p^J . 

and corresponds to a standard shower with vetoes applied such that only emissions with 
P± < P±h generated. The truncated shower is defined by 

Sf {qi,qh;P±iJ = ^Y{qi,qh)S~ {qo)+ / {qi,q)dVi-^ig{q, z)9{p^,^ - p±{q, z)) (3.10) 

Jqh 

xSf {zq,qh;p±J Sj ((1 - z)q;p^J . 

and corresponds to a standard vetoed parton shower line, constrained not to produce any 
flavour changing emissions, that is stopped once the truncated line has evolved down to 
the scale tjh. 

3.2 CKKW shower reorganisation 



In the POWHEG treatment reviewed in Sect. 3.1 a single hardest emission is separated such 



that it may be corrected with matrix elements. In the CKKW algorithm we aim to improve 
the parton shower with tree-level matrix elements for all parton multiplicities resolved at the 
merging scale A;j_^^g . We perform a reorganisation of the parton shower, analogous to the 
POWHEG reorganisation, splitting the shower into two parts: a hard shower describing 
emissions resolved above the merging scale; and another shower producing the rest of 
the shower emissions around this hard shower. The hard shower can then be generated 
according to the tree-level matrix elements as required by the CKKW algorithm. 

The result of this generalisation of the POWHEG reconstruction is a set of truncated 
and vetoed showers which fill in the radiation between the hard emissions defined by the 
hard shower history. In order to see how this works we first consider the next step up from 
the POWHEG case of a single hard emission, where we have exactly two hard emissions 
along the hard shower line, generated at scales qi and q2- One possible configuration of 
this hard shower line is given in Fig. |l[ As was done in formulating the POWHEG scheme, 
the full parton shower can be constructed around this hard shower line by constructing an 
equation analogous to Eq. (^), 

rqi 

S^^iqi) = / ST^ {qi, qi; k^,,s) .~(qi,Zi)SY ((1 - z,)qr, A:^,, J (3.11) 

ijk Jqo ijk ijk^ijk 

: z . u \ ^T>_ ( z ^ \ cV ( ^ z . u \ CfV ft -I ^ \z . u \ 



The superscript (2) on S denotes that this does not describe a general shower line, but 
the subset of shower lines with exactly two emissions above the merging scale. Eq. (|3Tll ) 



Ml 






Si, 


J go 


ijk 


Ml 










jk 







-12- 



Figure 1: An example of a hard shower hne configuration where two emissions are generated above 

k±MS ■ 



contains two truncated showers, one containing parton shower emissions with k± < k^j^^^ 
before the hard emission and the other containing emissions with k± < k±j^^g 

between the hard emissions at {qi^zi) and {q2,Z2)- The truncated and vetoed showers are 

respectively. 



given by Eq. (3^) and Eq. (|3l 

As in the POWHEG case, the splitting functions and Sudakov form factors for the 



truncated and vetoed showers in Eq. (3.11) do not match each other and are therefore not 
suitable for a standard Monte Carlo treatment. However, we can use the same manipula- 
tions as in the POWHEG formulation to split the Sudakov form factors into a product of 
a Sudakov form factor that matches the vetoed splitting functions and a remnant Sudakov 
form factor, as in Eq. (|3.5| ). The truncated showers contain only soft radiation so again 
we may set Zi ^ 1 in the remnant Sudakov form factor of Eq. (l3.7| ) with only sub-leading 
differences. The result of this is that the product of remnant Sudakov form factors for 
a particular truncated or vetoed line combine to give a remnant Sudakov factor. Rather 
than resulting in a single remnant Sudakov factor as in the POWHEG scheme, we now get 
a product of remnant Sudakov factors. The product of remnant Sudakov factors for the 
hard shower configuration of Fig. |l| is given by 



ijk 

X Af{z2q2; k_ 



-Af((l-zi)gi;A;x,«) 



(3.12) 



)A§iil-Z2)q2;ki_ 



MS Jl 



where the remnant Sudakov factor is given by 



A5(g; kj 



■MS I 



exp 



Jqa ~. . . 



{q,z)e{k±{q,z) - A;_L„s) 



(3.13) 



In the CKKW algorithm the hard shower is generated by choosing a jet multiplicity n 



as described in Sect. 2.1 and generating n parton momenta according to the appropriate 
jet cross section. A pseudo-shower history and corresponding shower variables are then as- 
signed by applying a clustering algorithm to the n parton momenta, until they are clustered 



®With the replacement p± ^ k±. The variable p± is the shower definition of transverse momentum and 
k± is the merging variable transverse momentum measure. 
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back to a leading order configuration. The shower reorganisation presented here results 
in a product of remnant Sudakov factors, with which these hard shower configurations 
should be reweighted. These remnant Sudakov factors can generally be found from the 
pseudo-shower history by applying the following prescription: 



each internal line from a branching at (qi, zi) to {q2, Z2) contributes a factor 



(3.14) 



• each external line from a branching at (q, z) contributes a factor 

Af(zg;/cx,«). (3.15) 

These remnant Sudakov factors match the Sudakov factors in Eq.( |2.1l| ) that we argued 
should be introduced in order to extend the CKKW procedure for transverse momentum 
showers to the angular-ordered shower. 

The parton shower for emissions below the cut is generated by producing truncated 
and vetoed showers around the hard shower according to the following prescription: 



• each internal line from a branching at (gi, zi) to (52, Z2) results in a truncated shower 

-5/ (2191,92; ^±Ms); (3-16) 

• each external line from a branching at (9, z) results in a vetoed shower 

5'/ (29;^±Ms) • (3-17) 

4. The algorithm 



In order to implement the procedure described in Sect. 3.2 we employ the strategy of 
Ref. [18], where the hardest emissions, or set of hard emissions in this case, are interpreted 
as parton shower emissions. This approach leads to a straightforward implementation 
of the truncated showers, where a truncated shower, evolving between hard emissions at 
(gi, zi, (f)i) and (^2; Z2-, 4>2)-, is generated by initiating a standard parton shower at zi^i with 
vetoes allowing only non-flavour-changing emissions with k±_ < k±j^^g and stopping the 
truncated shower once it has evolved beyond q2, at which point the second hard emission is 
forced with splitting variables (^2) -22; 02)- This allows the full shower of truncated showers, 
hard emissions and vetoed showers to be generated as a single shower evolution from the 
leading-order configuration. This results in a substantial improvement over earlier CKKW 
implementations with angular-ordered showers [7,12], since now the colour structure in the 
event is plainly equivalent to that which the shower would have produced by default i.e. 
it respects colour coherence. 

In order to interpret the matrix element emissions as shower emissions, we require an 
exact mapping from the set of n external parton momentum and assigned pseudo-shower 
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history to a set of shower splitting variables, {q, z, (j)), describing each emission. Obtaining 
such a mapping equates to inverting the momentum reconstruction, which is performed 
at the end of the standard parton shower to translate the set of shower variables into the 



parton momenta, and is described in Sect. 4.2. Having such a mapping also provides the 
exact shower variables that are to used for the Sudakov and as reweighting. 
The full modified CKKW algorithm is described below. 

1. The jet multiplicity n is generated with probability 

where cross sections are evaluated at a fixed strong coupling asj^j^ ■ 

2. The n external parton momenta are generated according to dcr(A;_Lj,^g). 

3. Pairs of external parton momenta are clustered'' down to a leading-order configura- 
tion, assigning a pseudo-shower history. 

4. The inverse momentum reconstruction is applied to the external momenta and shower 
history such that a set of shower splitting variables (g, z, (p) are found, describing n — 2 
hard branchings. 

5. The configuration is reweighted to include the Sudakov form factors and running 
as- This corresponds to assigning the configuration a weight W and rejecting the 
configuration if < Ti? . The weight is constructed from the pseudo-shower history, 
according to the following prescription: 

• each hard emission at (g, z) contributes a running as factor 

"5 



"5a 



(4.2) 



each internal line between hard emissions at {qi,zi) to (q'2,-22) contributes a 
Sudakov factor 

A^(2;i^i;A;_L^,g)_ 



(4.3) 



• each external line from a hard emission at (g, z) contributes a Sudakov factor 

^f{zq;k^,,s)- (4-4) 
If the configuration is rejected^ return to step ||. 



^The clustering procedure is discussed in Sect. 4.4 



TZ refers to a random number, generated in the interval [0, 1]. 
®This reweighting procedure relies on the weight generated in this step satisfying W < 1. The fixed 
strong coupling used in the matrix elements a^^^^ can be chosen to be large enough that the as weight is 
always less than one. Individual Sudakov form factors are also guaranteed to be less than one while the 
ratio of Sudakov form factors contributed by intermediate lines must be less than one due to the angular 
ordering condition ZiQi > qi+\. 
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6. Parton shower lines are initiated from the leading-order configuration which are to 
be evolved according to the procedure: 



(a) If a hard emission exists at a lower scale on the shower line, then the shower is 
evolved as a truncated shower otherwise proceed with step (|6^) . The truncated 
showers evolve as the standard parton shower with vetoes allowing only non- 
flavour-changing-emissions with k±_ < fc_L^^^ . Each truncated emission generates 
a soft gluon which should be evolved according to step (|6cD. 

(b) Once the scale of the next hard emission is reached, the hard emission is forced 
creating two further shower lines, each of which should be evolved according to 
step (|6al). 

(c) Vetoed showers evolve all external shower lines down to the hadronization scale, 
with vetoes allowing only emissions with k± < k^^j^. 

The above scheme is adapted for the highest multiplicity channel, where n = A'", by the 
replacement k±j^jg — > k_\_^ in the shower vetoes and Sudakov form factors. 

4.1 Shower kinematics 

In order to implement the inverse momentum reconstruction and parton shower vetoes, 
the kinematics of the Herwig-|— |- final-state parton shower must be understood. The full 
details of this are given in Ref. [2] and we present a review here for completeness and to 
introduce our notation. 

Each external parton from the hard sub-process with momentum pj is interpreted 
as a progenitor for a parton shower jet. Each parton shower jet evolves from the initial 
scale qi = y/s down to the hadronization scale qq undergoing a series of branchings, each 
described by the shower variables {q^z^cj)). 

Once all of the shower lines have evolved down to the hadronization scale, the shower 
evolution is stopped and the momentum of all external and intermediate partons are re- 
constructed from the shower variables. This is done in the centre-of-mass frame via the 
Sudakov decomposition, where the momentum of a parton in the shower is written 

qi = oiipj + f5inj + qu, (4.5) 

where the reference vector nj is taken to be a light-like vector with three-momentum equal 
to that of the colour partner to the jet progenitor and q± is the component of momentum 
transverse to both pj and nj. The momentum fraction z of each branching is defined by 

^ = — , (4.6) 



the scale of the emission q is defined in Eq. (2^). The relative transverse momentum of 
the branching is defined by 

P±i = q±i - zq^-j. (4.7) 
The p±i vector is written in terms of the azimuthal angle (j) 

P± = {\P±\ coscj), \p±\ sin(/),0;0) . (4.8) 
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The Sudakov variables Oj, Pi and q±i are calculated from the shower variables recursively 
for all partons in the shower jet and the momentum are constructed according to Eq. ([4.5|). 

The momentum reconstruction procedure results in the progenitor momenta qj gaining 
some off-shell momentum and leads to the loss of global momentum conservation. Lon- 
gitudinal boosts are applied to each shower jet to restore momentum conservation while 
disturbing the internal structure of each jet as little as possible. These reshuffling boosts 
are defined for each jet by the transformation 

V«1 + ij) ^ (^^^■^; y^^5+^) • (4-9) 

Momentum conservation is ensured by requiring that the rescaling parameter k satisfies 



E \/^'Pj + 1' = ^'- (4-10) 
J 

4.2 Inverse momentum reconstruction 

In the CKKW algorithm a clustering procedure is applied to a configuration of external 
parton momenta, defining a pseudo-shower history for producing that configuration as a 
set of branchings from a leading-order parton configuration, which in this case would be 
qq. We aim to interpret this branching history as a set of shower branchings such that 
when this set of branchings are forced, the momentum reconstruction will reconstruct the 
original parton momenta. In order to map the shower history to a set of shower variables 
the momentum reconstruction procedure must be inverted; this requires two steps. First, 
the boost applied to each shower jet in order to conserve global momentum must be found 
and its inverse applied to the momenta of the shower jet. Second, the resulting momenta 



are decomposed into the shower variables according to Eq. (4.5). 



The momenta of the set of progenitors defined by the branching q'j correspond to the 



momenta on the right hand side of Eq. (4.9). The original on-shell progenitors pj are 
related to q'j by 

where mj is the on-shell mass of the jet progenitor. The set of on-shell progenitors respect 
global momentum conservation therefore we can find the boost parameter k by solving 



J 



Once k is found, the reference vector pj is given by Eq. ( [4.11 ); similarly nj is given by 




where J refers to the colour partner jet of the shower jet J. In the reconstruction proce- 
dure, the Sudakov parameters of the progenitor partons are set to oq = 1 and q_\_Q = 0. 
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Furthermore, Eq. ( [4 .51 ) implies that 



/3o = (4.14) 
2p ■ n 



Since q'j and qj are related by a boost, we also have = q'J. The momentum of the 
reconstructed progenitors qj can then be constructed according to Eq. ([4.5|). This defines 
the reshuffling boost as in Eq. ( [4.91 ). The boosts for all shower jets can then be calculated, 
inverted and applied to all momenta in each jet. The momentum can then be decomposed 
into Sudakov parameters and the shower variables (g, z, (p) for each branching calculated 
from Eqs. (^I^-pl 



4.3 Shower vetoes 

The vetoes that are applied to the truncated and vetoed showers and the cuts applied to 
the remnant Sudakov form factors require a mapping between the shower variables, {q, z) 
and the merging scale transverse momentum measure k±. The merging variable, for an 
emission ij ij, is defined in some jet measure according to, k± = ^/yijs. We have 
implemented the merging algorithm with the Durham [23] and LUCLUS [24] jet measures, 
defined by 

2min (Ef,E^) 

ydur,, = ^ (1 - COS %) , (4.15) 



2{E,E,) 



2 



= , ' 2 (1 - COS 9ij) . (4.16) 

S {tji + £jj) 

In order to implement these vetoes a mapping between the shower variables {q, z, (p) 
and Hij, in the chosen jet measure, must be found. The Herwig++ shower produces off- 
shell intermediate states and therefore a set of boosts must be applied to each shower line 
in order to ensure momentum conservation. Since the boosts depend on the full shower 
history, an exact mapping between the shower variables and the merging variable cannot 
be found. We use a mapping that is exact for a single shower emission and should give a 
good approximation for larger numbers of emissions. For clarity in the following, we treat 
partons as massless while in our implementation parton masses are retained. 

In the case of e^e~ — > hadrons we have two shower progenitor partons q{pa) and q{pb)- 
The momenta of these progenitors, in the centre-of-mass frame, can be written, 

Pa,b = ^(0,0,±l;l). (4.17) 

We now consider a single gluon emission along the quark line such that q{pa) — *■ <l{<li)g{<l2)- 
This is the first emission so we have ai = z and 02 = 1 — z. The transverse momentum of 
the emitted partons are given by q±i2 — ^P±- The emitted partons are considered to be 
external partons and therefore their momenta should be set on-shell. The (5 parameters in 
Eq. ( |4.5| ) are therefore given by 

A = ^, (4.18) 

ZS 
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P2 



(4.19) 



{l-z)s- 

The vetoes should correspond to vetoes on the reshuffled momenta that have had the boosts, 
defined in Eq. (p^), applied to them. We should therefore solve Eq. ([4.10| ) and calculate 
these boosts before applying Eqs.( 05| , 4.16| ). The reconstructed progenitor momenta are 
given by qa = qi + Q2 foi' the quark jet and qb = Pb for the anti-quark jet. Inserting these 
momenta into Eq. ( 4.10| ) yields the solution 



k = l- y (4.20) 

szyV — z) 

for the boost parameter. The reshuffling boost for the quark line is then defined by 
Eq. (|4.9|). It follows that the three- vector of the shuffled quark progenitor should be 
given by 



(4.21) 

The expression in Eq.( 4.2l| ) is exactly the same as qa and therefore the boost to be applied 
to the quark jet is the unit matrix. The shuffled momenta for the emitted partons have 



now been constructed and we can apply Eqs. (4.15, 4.16) to give expressions for the jet 
measures used to define the merging scale. These give 

2 



2/du 



mm 



sz 



^ + 



s{\-z) 



1 



COS ( 



yiuc 



{v\ + (1 



p^s + s^z(\ — z) 



(1 -cos6') 



where 



cos I 



1 



2p2 s 



(4.22) 



(4.23) 



(4.24) 



{pI + (1 - z)2s) {pI + ^2^) • 

These mappings allow a transverse momentum measure, k± = ^Jys, to be calculated in the 
merging variable for each parton shower emission. Parton shower vetoes and Sudakov cuts 
can then be applied by comparing this measure to the merging scale. 

4.4 Clustering scheme 



The parton shower decomposition presented in Sect. 3.2 relied on our ability to interpret the 
series of hard branchings, defined by the matrix element momenta and assigned pseudo- 
shower history, as a parton shower. The inverse momentum reconstruction procedure 
ensures that, given an assigned pseudo-shower history, a set of parton shower emissions are 
found that will exactly reproduce the matrix element momenta. 

The assignment of a pseudo-shower history is dependent on the clustering scheme used 
and this should be carefully chosen to ensure that the assigned history is similar to that 



which the parton shower would have generated. In particular. Sect. assumes that the 
assigned history is angular-ordered. We therefore assign histories that obey the angular- 
ordering condition 

qiZi > qi+i, (4.25) 
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for all emissions along all shower lines. 

The inverse momentum reconstruction allows us to find the shower variables of all 
branchings in a particular pseudo-shower history. We can therefore determine whether a 
history is angular-ordered by following all shower lines outwards from the hard process and 



explicitly checking that all of the branchings satisfy Eq. (4.25). 

We employ a clustering procedure that creates all possible angular-ordered pseudo- 
shower histories and chooses the one that the shower was most likely to produce. This is 
achieved by the procedure: 

1. all possible shower histories are created by clustering all pairs of partons whose 
flavours correspond to allowed branchings; 

2. non-angular-ordered histories are discarded; 

3. the angular-ordered history for which the scalar sum of the transverse momentum of 
its branchings is smallest is chosen; 

4. if no angular-ordered histories could be found, the unordered history with smallest 
transverse momentum is chosen. 

The procedure of always choosing the shower history with smallest transverse momentum 
is taken rather than a probabilistic choice so that unnatural clusterings are not assigned, 
for example in a three jet event a gluon collinear to the quark line could be clustered to 
the anti-quark. 



5. Results 

In this section we present the results of the implementation of the modified CKKW al- 
gorithm for the process e~^e~ hadrons at a centre-of-mass energy of 91.2 GeV at both 
parton and hadron level. The parton level results provide a test of the algorithm's ability to 
provide a smooth merging between the matrix element and parton shower regions of phase 
space, showing features that may be hidden by the addition of a hadronization model. The 
hadron level results provide the ultimate test of the algorithm's ability to describe data 
and in particular are sensitive to the parton colour structure assignment which we expect 
the modified algorithm to improve with respect to traditional CKKW methods. 

A key test of the merging algorithm is its insensitivity to changes in the merging scale 
and merging variable. The algorithm was implemented with two merging variables: the 
Durham and LUCLUS jet resolution variables. For each merging variable, merging scales 
of yJ^^g = 5 X 10~^, Uj^jg = 10~^ and y^^^ = 5 x 10~^ were used. Samples of events with all 
partons separated hy y > y^j^ were generated using MadGraph/MadEvent [25,26] for the 
process with up to five partons in the final state^'^. 

^^The maximum multiplicity for each merging scale was decided according to the phase space available 
in the matrix element region. 
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Figure 2: Parton level distributions of the scale at which three jets are resolved in the Durham jet 
measure for e+e^ hadrons at y/s = 91.2 GcV. The red line shows the CKKW distribution with 
maximum multiplicity set to three, the blue line shows the Herwig++ parton shower distribution 
with a matrix element correction applied and the cyan line shows the HerwigH — h parton shower 
distribution without the matrix element correction. Plots (a)-(c) show the CKKW distributions 
with merging scales set to y^^^^ = 5 x 10^'^, yj^^,, — 10"^ and y^^g = 5 x lO"-^ in the Durham jet 
measure. Plot (d) shows the CKKW distributions at all of the merging scale choices on the same 
plot. 



5.1 Parton level results 

We present the distributions of the merging variable itself since these should be the most 
sensitive to problems with the merging procedure. In order to provide a direct comparison 
to Ref. [13], we first present a systematic look at the algorithm with the maximum multi- 
plicity set to three, so that the matrix element region is responsible for, at most, a single 
hard emission. 

Figure |2| shows distributions of the scale at which three jets are resolved in the Durham 
jet measure for the three chosen merging scales with the Durham jet measure as the merging 
variable. Jet analyses were performed with the KtJet package [28]. Each of the merging 
scale choices exhibit a smooth transition between the two phase space regions, there also 
appears to be little dependence on the choice of merging scale. 

Figure |3| shows the same distributions as Fig. || but with the truncated shower switched 
off. Switching off the truncated shower results in a radiation gap, meaning that emissions 
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that would be generated at scales greater than that of the hard emission, but with transverse 
momentum less than that of the hard emission, are never produced. This radiation gap 
corresponds to a deficit in the amount of soft wide angle emissions produced from the 
three jet samples. This problem will only become more serious as higher multiplicity 
contributions are included and underlines the importance of the truncated shower in the 
merging algorithm. 

Figure ^ shows the same distributions as Fig. |^ but with the highest multiplicity 
treatment switched off. The result of switching off the highest multiplicity treatment is 
that a maximum of three emissions may be generated in the matrix element region. This 
violates the all-orders-in-a^ resummation of the parton shower. The result is a suppression 
of the three-jet channel and a deficit in the hard radiation generated in the three jet channel. 
In Fig. ^ this presents itself as distributions that are peaked around the merging scale and 
have a large dependence on the choice of merging scale. 

Figures |5| and ^ show the distributions of the scale at which three jets are resolved, 
for the algorithm with maximum multiplicity set to up to five jets, with the merging 
algorithm defined in the Durham and LUCLUS jet measures respectively. As in Fig. |2|, all 
distributions appear to be smooth around the merging scale and to be relatively insensitive 
to the choice of merging scale. 
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Figure 4: The same distributions of Fig. || but witli the highest multiplicity treatment switched 
off in the CKKW treatment. 



Since we are now including higher multiplicity channels in our merging algorithm we 
check the distributions of scales at which higher numbers of jets are resolved. This is done 
in Fig. ^ for the resolution of four and five jets in the Durham and LUCLUS jet measures. 
The merging in these distributions is well behaved. 

These distributions demonstrate a degree of insensitivity to the choice of merging 
scale, which has been varied over an order of magnitude, however there is still some residual 
dependence on this choice. While the parton shower and merged matrix element treatments 
formally have the same large logarithm behavior, there are differences between the two. 
The degree of these differences will directly influence the amount of residual dependence 
on the merging scale that is observed. In changing the merging scale we are changing the 
volume of the matrix element phase space region and therefore changing the proportion of 
parton emissions that are corrected by exact matrix elements. 

It can be seen in Fig. ^ and even more so in Fig. ^ that when the merging scale is set 
to its upper value of y^,^^ = 5 x 10^^, the volume of phase space into which the matrix 
elements are allowed to emit is small. This results in the CKKW distributions being very 
close to those of the Herwig++ parton shower with no matrix element correction. Lowering 
the merging scale increases the volume of the phase space region described by exact matrix 
elements and moves the distributions farther from those of the bare parton shower. Table 
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Figure 5: Distributions of the scale at which three jets are resolved in the Durham jet measure. 
The red line shows the distributions for the CKKW treatment with all multiplicity channels (up to 
a maximum of five jets) included at a set of merging scale choices in the Durham jet measure. 



|l] gives the cross sections for the CKKW treatment at different choices of the merging scale 
and exhibits variation at the 5% level. 



Vms 


Durham cross section / nb 


LUCLUS cross section / nb 


5 X 10-2 


38.2 


38.6 


10-2 


36.5 


37.1 


5 X 10-3 


35.7 


35.9 



Table 1: Table of cross sections of the process e+e hadrons for different choices of the merging 
scale in the Durham and LUCLUS jet measures. 



5.2 Hadron level results 

We present a comparison of the Herwig++ CKKW implementation with hadronization 
switched on to LEP data for a variety of event shapes. It is standard practice to tune the 
free parameters of an event generator to LEP data and this has been done with the default 
Herwig++ parton shower with matrix element corrections. Since the CKKW merging 
algorithm significantly changes the parton shower component of the event generator and in 
order to provide a fair comparison with default Herwig++, a new tune was performed on 
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Figure 6: Distributions of the scale at which three jets are resolved in the LUCLUS jet measure. 
The red line shows the distributions for the CKKW treatment with all multiplicity channels (up to 
a maximum of five jets) included at a set of merging scale choices in the LUCLUS jet measure. 



the free parameters for Herwig++ with the CKKW algorithm. This tune was performed 



with the merging scale set to y^^^g = 10~ in the Durham jet measure. The most notable 
change to the default Herwig++ parameter set was a change in the value of the strong 
coupling at Z mass from a value of as{M2,) = 0.128 for Herwig++ with matrix elements 
corrections down to as{Mz) = 0.120 for Herwig++ with CKKW. This change brings the 
parameter closer to its measured value and is indicative of an improved treatment of hard 
radiation. 

Figures |-10 show distributions of a range of event shape, jet resolution and four-jet 
observables in comparison to LEP data. The parton level analysis shows that the merging 
scale choice of y^^g = 5 x 10~^ leaves only a very small region of phase space that is corrected 
by the matrix elements. This very high scale choice will therefore not give the improvement 
expected in introducing the merging algorithm, we therefore omit this merging scale choice 
from the hadron level analysis. In each of the figures the red band shows the variation in 
distributions over the four merging scale choices of y^^g = 10^^ and yj^^g = 5 x 10^^ in the 
Durham and LUCLUS jet measures. 

The CKKW distributions (red band) in Figs. |^-|lO| all demonstrate improved descrip- 
tions of the data in comparison to the default Herwig++ parton shower with matrix element 
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corrections. In particular the tails of the distributions in Fig. ^, corresponding to hard emis- 
sions, and the jet resolution distributions of Fig. |9| with four and five jets are significantly 
improved as would be expected given the aims of the merging algorithm. The four-jet angle 
distributions of Fig. |l^ are also all improved, with the exception of the 034 angle, which 
was already well described by the default Herwig-|— |- parton shower. The Onr distribution 
provides the most notable improvement in its description of the data in comparison to the 
default Herwig-|— I- parton shower. 

The width of the red band on the distributions shows that there is some residual 
dependence on the merging scale however it does not appear to be too serious and is at a 
similar level to that observed at parton level. This shows that the problems with colour 
structure, that appear in the standard CKKW algorithm, are not present here and that the 
truncated shower is working as intended. It should be noted that a fixed set of Herwig-|--|- 
shower and hadronization parameters was used for each of the four merging scale choices; 
the variation would be reduced further if a tune of the parameters was performed for each 
merging scale choice. 

The per degree of freedom values for the distributions in Figs. 0-10 are given in 



Table | for the merging scale choice of y^^-g = 10 ^ in the Durham jet measure, which was 
used in the tune. The CKKW values are lower than those of the default Herwig-|--|- shower 



-26- 



Oblateness 
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Sphericity Planarity 
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Figure 8: Distributions of the event shape variables thrust, oblateness, sphericity and planarity 
for e'^e" hadrons at a centre-of-mass energy of y/s = 91.2 GeV in comparison to LEP data 
(black) [32]. The red band gives the variation of the distributions of the CKKW implementation 
with merging scales choices of y^^j^ — 10~^ and 2/^,^3 = 5 x 10^"^ in the Durham and LUCLUS jet 
measures. The blue histogram gives the distributions of the default Herwig++ parton shower with 
matrix element corrections. The lower panel shows the ratio of the difference between simulation 
and data to the data in comparison to the error bounds of the data (yellow region) . 

in all cases except for the 034 angle, where the default implementation already gave an 
accurate description, and in many cases the CKKW values are significantly lower. 

6. Conclusions 

A modified version of the CKKW algorithm has been implemented in Herwig++ for the 
process e~^e~ — > hadrons. The modified algorithm uses truncated showers in order to pro- 
vide smooth merging between the Herwig-|— |- angular-ordered parton shower and a set of 
transverse-momentum-ordered emissions defined by inverting the Herwig-|— |- momentum 
reconstruction procedure on a samples of parton momenta generated according to exact 
tree-level matrix elements. 

The truncated shower was found to result in a smooth merging between parton shower 
and matrix element regions of phase space with parton level distributions appearing free 
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Figure 9: Distributions of the scale at which three, four and five jets are resolved in the Durham 
jet measure for e+e" hadrons at a centre-of-mass energy of y/s = 91.2 GeV in comparison to 
LEP data [33]. The colours of the lines are the same as those in Fig. ||. 

of discontinuities around the merging scale and relatively insensitive to changes in the 
merging scale. 

A full tune of the Herwig++ free parameters was performed for the CKKW implemen- 
tation with a merging scale of y^^^ = 10^^ in the Durham jet measure. This was found 
to give a good description of LEP data, demonstrating a significant improvement over the 
results from the default Herwig-|--|- parton shower with matrix element corrections applied. 

The results show a comparable level of merging scale dependence and agreement with 
LEP data to that found in Ref. [20], in which a similar CKKW merging approach was 
performed with a transverse-momentum-ordered dipole shower. 
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